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We investigate the isotope effect on the London penetration depth of a superconductor which 
measures ns/m*, the ratio of superfiuid density to effective mass. We use a simplified model of 
electrons weakly coupled to a single phonon frequency we, but assume that the energy gap A does 
not have any isotope effect. Nevertheless we find an isotope effect for ns/m* which is significant 
if A is sufficiently large that it becomes comparable to we, a regime of interest to high T c cuprate 
superconductors and possibly other families of unconventional superconductors with relatively high 
T c . Our model is too simple to describe the cuprates and it gives the wrong sign of the isotope 
effect when compared with experiment, but it is a proof of principle that the isotope effect exists for 
ns/m* in materials where the pairing gap and T c is not of phonon origin and has no isotope effect. 

PACS numbers: 74.25.Ha, 74.25.fc, 74.25.Kc 



I. INTRODUCTION 

While there is a general consensus that strong correla- 
tion governs the basic physics of high T c cuprates^, the 
role of electron-phonon interaction in determining T c is 
still under debate. The isotope effect is often viewed as an 
useful tool which can provide information on this impor- 
tant issue. There is extensive isotope effect data on hole 
doped cuprates and the following picture has emerged. 
There exists significant isotope effect on the transition 
temperature for underdoped cuprates, but none for over- 
doped ones. On the other hand, for all doping sub- 
stantial isotope effect on the London penetration depth 
X a b has been observed 2 ^. Recall that A~ b 2 oc ns/m* 
is a direct measure of the ratio between the superfiuid 
density ns and the carrier mass to*. The unusual iso- 
tope effect on T c can be understood qualitatively by the 
following picture^. For underdoped samples the tran- 
sition temperature is controlled by the phase stiffness 
A'5 = h 2 ns/4m*—. Hence, the isotope effect on T c may 
be simply inherited from the isotope effect on ns /to* . On 
the other hand, in overdoped samples T c is controlled by 
the pairing gap which has no isotope effect if phonons do 
not contribute significantly to its origin. Thus the isotope 
effect of T c can be qualitatively understood provided we 
accept the isotope effect on ns/m* . Then the puzzle is 
shifted to the origin of the isotope effect on ns /m* . Up to 
now there has been very little discussion in the literature 
on the isotope effect on ns/m* in superconductors. In 
this paper we take the first step to address this important 
issue. 

In 2 ^ the authors suggested that the isotope effect 
on ns/m* is due to an isotope effect on the effec- 
tive mass to*. However, the effective mass is given by 
771* /m = 1 + A where A is the dimensionless electron- 
phonon coupling. Usually A is not considered to have 
an isotope effect^. This is because A = g 2 N(Q)ujo and 
g 2 w (/ 2 ) /ujd, where / is the coupling of the electron 
density to lattice deformation, which has no isotope de- 
pendence. This is supported by the direct experimen- 



tal observation of Iwasawa et al£ Using high-resolution 
laser ARPES measurements of electronic dispersion in 
Bi2Sr2CaCu208+i (Bi2212), these authors study the in- 
fluence of isotope oxygen substitution on the boson cou- 
pling "kink" in the electronic dispersion. The experiment 
clearly reveals an isotopic dependence of the kink energy, 
while it observes no change in the effective electron mass. 
On the other hand this experiment provides evidence for 
the coupling of phonon mode£ to electrons, which leads 
us to propose alternative explanation for the isotope ef- 
fect on ns- 

Our model disregards any effects of strong correlation 
and considers a clean superconductor with s or d-wave 
singlet pairing. We suppose that there exists one phonon 
mode weakly coupled to the electron system. However, 
the superconducting order parameter is not related to the 
phonon mode, but rather induced by some other mecha- 
nism. Therefore, we ignore the isotope effect on the gap 
A and T c , even though these effects can easily be added 
to lowest order in A. We assume the standard adiabatic 
approximation for electron-phonon coupling^ and neglect 
Coulomb interaction. The only quantity which depends 
on the isotope mass is we. 

First, we present simple qualitative arguments that ex- 
plain the existence of the isotope effect on 775 due to 
phonons using the sum rule which relates total spectral 
weight (SW) to the electron density and (bare) electron 
mass, 

= - / duja (w), (1) 

TO 7T J Q 

where cr'(w) is the dissipative part of the conductivity. 
The onset of superconductivity does not change the to- 
tal SW, but is redistributed^. Some part of the SW now 
goes to the 5-function response of the condensate, whose 
weight determines ns/m*, while the remaining goes to 
frequencies w > 2 A. This is visualized in the FIG. [U 
It displays the dissipative part of conductivity for nor- 
mal metal, er'(w), that consists of a narrow Drude peak 
(since our sample is clean) along with a contribution from 
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FIG. 1. Qualitative form of the conductivity for a normal sys- 
tem. With the onset of superconductivity, all spectral weight 
to the left of the corresponding red line goes to the condensate 
response. 

phonons at some characteristic frequency w^. Qualita- 
tively with the onset of superconductivity some SW from 
frequencies ui < A goes to the <5-peak and the SW from 
frequencies u> > A goes to the optical conductivity of 
superconductor. Thus, if we have A <C uje (the case in 
weak coupling BCS theory), variation of uje would have 
no effect on SW that goes to the condensate response. 
This may be the reason why this problem has not re- 
ceived much attention in the literature up to now. If A 
and uje have the same order of magnitude, ns/m* re- 
sponds to the change of uje and isotope effect for ns/m* 
is present. 

In the following sections we present calculations for the 
isotope effect for ns/m* based on the sum rule. In Sec- 
tion |n] we introduce our toy model, then calculate con- 
ductivity and supcrfiuid density. In Scction lTni we discuss 
results and various approximations that were made. Fi- 
nally, in the Appendix we present detailed calculations 
of the response kernel and conductivity. 



II. CALCULATION FOR s AND d-WAVE 
SUPERCONDUCTORS 

In this section we present our calculation of the isotope 
effect on the superfluid density ns- We note, that both 
the experiment on the London penetration depth and the 
theory described in Eq. ([T]) give the combination ns/m* . 
From hereon we simply define ns as a product of ns/m* 
and m* which is taken to be constant. We calculate the 
dissipative part of the conductivity due to phonons. After 
this, having the expression for the conductivity via the 
sum rule, we get access to the superfluid density. 

The theoretical model is a superconductor with 
electron-phonon interaction, described by the standard 
Frolich Hamiltonian. The superconducting order param- 
eter (singlet pairing with s or d wave gap symmetry) is 
assumed to be not related to phonon mode under consid- 
eration. One may think about this pairing as induced by 
other phonon modes, strong correlation or some differ- 
ent mechanism. The resulting Hamiltonian can be writ- 



ten as: 

H = Hbcs + Hph + H c . p h. (2) 

Hbcs is the BCS Hamiltonian, with the pairing 
in Cooper channel treated within mean-field theory 
(Coulomb interaction is neglected): 

#bcs ^J2^p t ^p 
p 

-^E^( A p^+ A p t -)*P' ( 3 ) 
p 

where £ p is the Bloch energy measured relative to the 
Fermi energy (we consider only one electron band that 
is coupled to one phonon mode), t 3 is the Pauli matrix. 
Electron field operators here are written as a vector in 
the Nambu space 



with c Pj ^, cp"^. being standard creation and annihilation 
operators for up (down) spin electrons. The Hamiltonian 
that describes the phonon mode and electron-phonon in- 
teraction is written as 

H p h + #e-ph = w ( k ) & k & k 
k 

+ E s(p, p'xv-p + K-p'K' t ^p- ( 5 ) 

p,p' 

Here w(k) is the (bare) phonon dispersion and <?(p,p') 
is the electron-phonon coupling. In what follows we use 
the Einstein approximation for the dispersion of phonons 

w(k) = const = uje, (6) 

and assume that our characteristic frequency is on 
order of the superconducting gap, i.e. uie ~ A (note 
that for BCS superconductivity u>e ^> A). Moreover, 
we limit ourselves to leading order corrections in the 
electron-phonon coupling that is assumed to be weak. 
Thus we can neglect not only renormalization of the 
electron-phonon coupling (that is small in the parameter 
uie/ Ep^-), but renormalization of the phonon dispersion 
as well. 



A. Green's function and self-energy 

Before we start our calculations of conductivity we 
need to consider the electron Green's function. Follow- 
ing standard procedure, we define the Green's function 
of electrons G p (t) (which is a matrix in Nambu space) 
and phonons Dq{t) in imaginary time as 

G p (r) = -(T T (* p (r)*+(0))> I (7) 
D q (r) = -(r T (6+(r)6p(0))>. (8) 
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FIG. 2. Diagrams corresponding to the self-energy (a) correc- 
tions and conductivity (b). 



Fourier transform of bare Green's functions (which are 
denoted by index (0)) can be written in the forroii: 

(G° {iujm))- 1 = lU) m - £ q T 3 - An, (9) 

Dl(iu v ) = - . 1 (10) 

4 — w(q) iuj u + uj{q) 

where ioj m = (2n + 1)ttT and iuj v = 2vkT are the Mat- 
subara frequencies for fermions and bosons. 

In order to find the Green's function of electrons in 
the presence of the electron-phonon interaction we note 
that the superconducting order parameter is not related 
to the phonon mode under consideration. Moreover, the 
electron-phonon coupling is small and we limit ourselves 
to the leading order correction. Therefore, we can use 
perturbation theory to get access to the electron Green's 
function and there is no need in doing the full Eliashbcrg 
strong coupling theoryi^. Renormalization of the Green's 
function is given by the self-energy S ph due to electron- 
phonon interaction: 

(GpHr 1 = (Gp^Hr 1 - s ph(p^) 

= Q r (uj)-£ p t 3 -A r (uj)t u (11) 

where w(w) and A(ui) are the renormalized uj and gap A. 

It suffices to calculate the self-energy to the leading 
order in electron-phonon coupling. The expression for 
the self-energy in the Matsubara diagram technique can 
be read from the diagram depicted in FIG. [2] (a): 

E ph (p,io;„) = f (dp')\g(p,p')\ 2 

D° p _ p ,(iuj n - iuj m )T 3 G a p ,(iuj m )T 3 . (12) 

Here the Pauli matrix T3 (note that the self-energy is 
a matrix in the Nambu space) comes from the interac- 
tion vertex. Since we work to leading order of perturba- 
tion theory, we use the bare electron and phonon Green's 
functions. We proceed with analytical continuation (the 
detailed procedure is described i n 12 i 13 ). Sums over posi- 
tive and negative m are represented as contour integrals. 
Afterwards, employing the spectral representation for the 
electron and phonon Green's functions (see booki 3 .), one 



has: 



sfh A ( P ^) 



d£ p , Im[T 3 G°?(z)T 3 ] 



tanh -Sf — coth 7^ tanh ^ + coth 



X — UJ — fi =p 16 



x — ui + n =p id 



(13) 



where a 2 (il)F(tt) is the spectral function of electron- 
phonon interaction (Eliashbcrg function) that consists 
of the product of effective electron-phonon coupling and 
phonon density of states. It is expressed through the 
electron-phonon interaction vertex as 



a 2 (n)F(n) 



Ifs d2 P 



(14) 



where v' F is the Fermi velocity and integrations are over 
the Fermi surface. Finally, we integrate over the loop 
momentum^: 



J dt p , lm[T 3 G° p R (x)r 3 } 



-i7rsignx 



y/x 2 - A 2 ' 



(15) 



and put temperature T = 0. Expanding S pl J (p, uj) in 
components we obtain for renormalized uj and A, that 
were defined earlier in Eq. (|lll) . 



r.R,A 



(uj) = UJ + 



dxdSla 2 F{Q)Rc 



x sign x 



V 1 x 1 



A 2 



0(-x) 



9{x) 



x — uj — Q =p i5 z — uj + Q id 



(16) 



A ra (uj) = A + dxdna 2 F(n)Rc 

J — OO 

6(-x) 6{x 



A sign x 



Q^f iS z — uj + fl iS 



(17) 



It is convenient to introduce a dimcnsionless coupling 
constant A defined to be responsible for the electron mass 
renormalization due to phonons: 



m* = m(l + A + 0(A)), 



(18) 



where m is bare and m* is renormalized (measured in 
experiments, e.g. in ARPES^) mass. The dimension- 
less coupling is expressed through the Eliashberg func- 
tion dH]) as: 



A = 2 / —a 2 (uj)F(uj). 
J u 



(19) 



We work in the Einstein approximation, where the dis- 
persion of phonon mode does not depend on the mo- 
mentum, uj p = ujfi. Using this, one can immediately in- 
fer from Eq. (jl"4")l that Eliashbcrg function a 2 (fl)F(Q) oc 
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6 (CI — uje)- The constant of proportionality can be read 
from Eq. (flU| , giving us the expression for the Eliashberg 
function that will be used in the remainder of this paper: 



a 2 (Cl)F(Cl) = ^8(Cl-w E ). 



(20) 



In what follows we will make use of the imaginary part of 
the renormalized frequency and gap, that can be easily 
inferred from Eqs. (fT%|) - (fT7)t ; 

-ft/// A m , , g(M -UE- A)(\uj\ -uj e ) (0 , 
uj (uj) = ttuj e \ 7= ; , (21) 



2y/(\ U \ -uje? 



A R " 



__,.^g(M-^-A)A S igna; 



2 A /(M-o; r ,)2-A2 

whereas, for advanced functions we have 

= -w*"(w), A a "(uj) = -A r "(uj). (23) 

Not surpisingly, the imaginary part of the self-energy is 
present only for uj > uje + A, i.e. where the real exita- 
tion can be created. The square root singularity that is 
present at this threshold would be smeared for a more 
realistic phonon spectrum. 

B. Conductivity 

Having the electron Green's functions at hand, we pro- 
ceed to the calculations of the conductivity. Following 
the standard approach^, in order to calculate the con- 
ductivity we consider the response kernel Q a p(k, ui) that 
relates current response to the vector potential: 

2 

TIG 

ja(k,w) = Q aP (k,u)A p (}i,u). (24) 

m 

The optical conductivity is expressed through Q as: 

ne 2 Q a p(k,uj) 



a a p(k,uj) 



m iuj 



(25) 



The response kernel Q a p(\t, uj) can be calculated within 
the Kubo linear response method as a sum of the diamag- 
netic contribution and a current-current correlator—: 



ne 



( k;W ) = S a0 -P R (k,oj), (26) 



where V R a(k, uj) is analytic continuation to the real fre- 
quency of the (Fourier transformed) current-current cor- 
relator 

V afj (v~v\r-T') = (r(j lQ (r,r)j 1/3 (r',r'))) . (27) 

The current operator ji(r) is defined as a paramagnetic 
part of the full current operator j (r) : 



j(r,r) = ^(V r 
2 m 



c+(r')c(r) 



-A(r)c+(r)c(r) = )x(v,t) - - A(r)c+(r)c(r). (28) 



Since we consider the superconductor in the London 
limit, we can put k = 0. Denoting Q xx (0,uj) = Q(oS) 
we proceed to the calculation of Q(oS). 

The current-current correlator V R a(k,uj) is given by 
the sum of all possible diagrams with 2 external vertices 
corresponding to current operators. If we neglect vertex 
corrections (see discussion in Section [TTTJ) , in the leading 
order in the electron-phonon coupling we need to consider 
only the simplest diagram, FIG.[5](b). It is convenient to 
represent the diamagnetic contribution to Q(uj) through 
Green's functions of a normal meta l 11 ' 15 . Combining 
the diamagnetic term with the current-current correla- 
tor given by the diagram in FIG. [2] (b) we have: 



Q(iu) n ) = - tr 



T^y] / d£ p (G p (iuj n + iuj m )G p (iuj m ) 
G£=°(iuj n + iu m )G$=°(iu m ))] . (29) 



Starting from this expression we perform analytic con- 
tinuation, and integrate over the loop momentum. An- 
alytical continuation is done in the standard way and is 
presented in details in the literatur o 13 ' 15 and in the Ap- 
pendix. After lengthy but straightforward calculations 
we get the following expression for the imaginary part of 
the response kernel Q(uj) at zero temperature: 



■ Im 

2 ,/A- 



dz 



-,AA 



+ 



(z + U), z) — 1 
' e A (z + uj)+e A (z) 
g RA (z + uj,z) + 1 ' 



e R (z + uj) - e A (z) 
Where e a (uj) with a = R, A is defined as: 



e a (z) = signz^(Cj a (z)) 2 - (A a (z)) 2 , 
and g a/3 are the structure factors introduced ir>i£: 

Q a (zi)Q p (z 2 ) + A a (z 1 )A' 3 (z 2 ) 



9 ap (zuz2) 



e a (z 1 )eP(z 2 ) 



(30) 



(31) 



(32) 



Using the expressions for the imaginary part of uj R ' A and 
A R ' A , Eqs. <J2TJ)-(J22J), we expand to leading order in the 
corrections to the Green's function. After some calcula- 
tions shown in the Appendix we arrive at the final an- 
swer, that can be expressed via complete elliptic integrals 
of the second kind, E(x) and K(x): 



ne 2 ttXujeuJ- 



E ^ - 



A(uj - uj e ) 



in 



m 



a'(uj) = 



where we used shorthand notation 

uj± = uj — uje ± 2A. 



K 



(33) 



(34) 
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FIG. 3. Correction to the dissipative part of conductivity 
cr'(w) due to phonons for normal metal (blue) and supercon- 
ductor with s or d wave symmetry (red and green curves corre- 
spondingly) for uje = 1.5A. Sigma normalized to ao = ^—X, 
Drude peak and (5-function response of condensate are not 
shown. 



In case when there is no superconductivity in our system, 
putting A = from (|33[) we reproduce the well-known 
result for the leading order correction to the conductivity 
of normal metal due to phonons^: 



°N,e-ph( w ) 



ne 2 ^itu!e(w — we) 



6(lu-uj e ). (35) 



ct'(w) has been computed earlier by Allenil. However the 
second term in Eq. (|33[) containing the elliptic function K 
is missing in Allen's formula. This difference is negligible 
when A <^ uje but gives a small noticeable correction 
when A > u)e- 



Gap renormalization A R (ui) — A = vanishes under av- 
eraging over angle, while the imaginary part of the renor- 
malization of u) reads: 



u) R "(u)) = ^9(\u)\-u) E ) 



M - u E \ 

V(M-w B ) 2 -A2/ 2 / FS 



7TU E X f\u)\-u) E 

—— 9(\u)\ - U) E ) K ' 



where k(x) is defined as 



, * 2 \xK(x 2 ) for x < 1 
w it I K{x- 2 ) for x > 1 



(38) 



(39) 



We repeat calculations of conductivity with A(u>) = 
A and the new expression for u) R "(uj). Correction to 
the conductivity to leading order in the electron-phonon 
coupling is 



a M = — ^3 



ne 2 2 / [-*\*\ zu r "(z + uj) 



dz ■ 



A|/|-a 



Vz 2 - A 2 



(40) 



FS 



Interchanging averaging over FS and integration over fre- 
quency we have two different expressions for loe > A and 
uje < A. When > A we obtain: 



a'(uj) = 



ne 2 ttujeX 



9{uj - uje) 



While for uje < A we have a more complicated expression 
involving the elliptic integral of the first kind F(sin$>; k): 



C. d-wave symmetry 

The previous calculations can be easily generalized for 
the d-wave symmetry of pairing. We use the simplest 
possible model with the cylindrical Fermi surface with 
an axis parallel to the c-axis of a crystal. The order 
parameter is supposed to have d x i_ y i structure. It can 
be written as 



A(k) = A/ with / = cos 26», 



(36) 



where 9 is the angle between momentum (which lies in 
the a — b plane) and a axis. 

First we recalculate our corrections to the electron 
Green's function. In contrary to the s-wave case, the 
gap A(k) is dependent on the direction of momentum. 
Therefore the previous answer for the s-wave case needs 
to be averaged over the Fermi surface. Averaging is de- 
fined as 



, ne 2 tkoeX 
a (uj) = — 0(b) - ui E ) 



m uj 
■o 



A-u 



{ \ z + U )\-u) E \J\^ 



(42) 



D. Isotope effect for superfluid density ns 



Using our explicit results for the conductivity, by 
means of the sum rule, we get access to the superfluid 
density. Initially the sum rule is formulated as 



ne 
m 



dui <j'{u)), 



(43) 



(X) FS = (d6) X(\cos2e\). 



(37) 



where n and m are mass and concentration of electrons, 
while cr'(ui) is the dissipative part of the conductivity. 
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FIG. 4. Correction to superfluid density for superconductor 
with s-wave (blue) or d-wave (red) symmetry. 
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FIG. 5. Isotope coefficient p defined in main text as a function 
of cue in our model. 



Let us first apply the sum rule to our system when 
there is no superconductivity, i.e. A = 0. We assume 
that the mean free path is big, thus scattering time due 
to impurities r satisfies tlue 1. Then, in addition to 
the very narrow Drude peak at w = 0, ctd(w), we have the 
small correction due to phonons, cr{ vrc ph (o;), that starts 
at wb (see Eq. (|35p). The calculation of the spectral 
weight of the Drude peak from the sum rule goes as fol- 
lows: 



this is usually the atomic mass of oxygen) was assumed: 



2 1 ~ a > i \ 



ne 
m 



X- 



dwcr N,o-ph(w) 



-(1 + 0(A 2 )), (44) 



where to* = m(l + A), see Eq. (|18]l. We see that the SW 
of CNe-ph( w ) docs not depend on and the sum rule 
reproduces the correct answer for the spectral weight of 
the Drude peak with the renormalized mass. 

With the onset of the superconductivity the picture 
does not change drastically. The dissipative part of the 
conductivity consists of the 5-function with the weight 
given by rise 2 /m* (response of the condensate) and the 
small contribution due to electron-phonon interaction 
c/(w) (calculated in Eq. ([33]) and Eqs. (gU)-(H3) for s 
and d-wave pairing). The sum rule implies for n<?: 



n s 



— ~ / auj a (w). 



(45) 



From Eqs. (|33|) and (j4lj) . the superfluid density ns can be 
easily computed. It depends onuj, thus giving nonzero 
isotope effect on superconducting density. Results of the 
numerical computations of (ns — n)/A as a function of 
wb/A are shown in FIG.IU 

Based on the dependence of ns on loe, an isotope co- 
efficient (3 can be easily calculated. It is defined as 2 - 



P = - 



(46) 



dlogAfo 2 d log uie 

where the usual relation between the phonon character- 
istic frequency and the atomic mass Mq (in experiment 



uje 



1 AM 

2 Mq ' 



(47) 



Results for (3 are shown in the FIG. [5J Notably, for both 
s and d-wave symmetry, the isotope coefficient has max- 
imum modulus when uje ~ A. While when loe ^> A, 
(3 vanishes in either case. Such behavior naturally fol- 
lows from our qualitative arguments, presented in the 
Introduction. However note that /? approaches zero from 
positive and negative values for the s-wave and d-wave 
pairing respectively. 



III. DISCUSSION 

Although our toy model proves that isotope effect on 
the superfluid density can exist due to phonons, it gives 
an effect of the opposite sign to the experimentally mea- 
sured. For example ir>2 for YPrBaCuO (3 = 0.38 for 
x = 0.3 and /3 = 0.71 for x = 0.4, while our model 
predicts the isotope coefficient of order of (3 ~ —0.1 for 
A~ 1. 

The opposite sign of isotope coefficient leads us to a 
discussion of the limitations of our model. We treated 
the electron-phonon interaction as being weak which is 
not really the case. As it can be inferred from the exper- 
iment^, the coupling constant is of order one. However, 
the qualitative features of our result should not change. 
The second approximation was the neglect of vertex cor- 
rections. The justification of this has been given mil. 
Since for normal metal ladder diagrams give minor ef- 
fects, we expect this to be the case for superconductors 
as well. 

The most serious limitation of our model with respect 
to high T c superconductors is the neglect of strong cor- 
relation effects. The minimum requirement would be to 
reproduce the fact that ns is proportional to the doping 
hole concentration x due to the proximity to the Mott 
transition, and not to the electron density 1 — x as in 
our toy model. Furthermore, from the point of view of 
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FIG. 6. Deformation of integration contour during analytical continuation. Branch cuts of Green's functions are shown in red 
and starts at lja and oja — iu) m - When we do analytical continuation from uj n — t —iuj, all four branch cuts lie on the real axis. 



the sum rule argument, it is known that strong corre- 
lation gives rise to an incoherent background at finite 
frequencies, the so-called mid-infrared peak, even before 
phonons are taken into account. This incoherent back- 
ground is missing in our toy model. What is needed is 
clearly a model which includes both Coulomb repulsion 
and coupling to phonons, such as the Hubbard-Holstein 
model. There has been considerable progress in this dif- 
ficult problem, as given in a recent review^. It may be 
possible to extend the present model to include correla- 
tion effect in the future. Even though our result does 
not agree with the experiment, the important message 
of the present paper is a proof of principle that isotope 
effect in ns /m* is possible even if the pairing is not due 
to phonons. 



Appendix A: Derivation of the response kernel and 
conductivity 

In this Appendix, followin g 11 ^ 5 we derive the response 
kernel Q(ui) and dissipative part of the optical conductiv- 
ity cr'(cj) due to the electron-phonon interaction. Start- 
ing from the expression for Q(iuj n ), Eq. we do an- 
alytic continuation to the real frequencies. Performing 
integration over £ p and setting temperature to zero, we 
arrive to the final expression for the response kernel. Af- 
terward, using corrections to the self-energy and expand- 
ing in the electron-phonon coupling we get the conduc- 
tivity. 

Matsubara sum in (|29p is represented as a contour in- 
tegral in the standard way: 



Finally the isotope effect on the supcrfhiid density dis- 
cussed in this paper should be applicable to conventional 
superconductors in the clean limit with minor extensions. 
To linear order in A, it is straightforward to include cou- 
pling to a distribution of phonon modes. Furthermore, 
we can easily include the isotope effect on the energy gap 
A in our consideration. In weak coupling BCS theory, A 
is proportional to uje and the ratio A/uje has no isotope 
dependence. Then our theory trivially predicts no iso- 
tope effect on ns/m*. Including /x* = ji/ \og(ep /uje) in 
the T c formula will introduce some isotope dependence 
in A/u)e- However, in most superconductors lue/A 3> 1 
and from FIG. [5] we see that the predicted effect is very 
small. Recently s-wave superconductors with relatively 
large T c and energy gap have been discovered. Examples 
are MgB 2 and doped fullerenes. While coupling to cer- 
tain high frequency phonons may be responsible for the 
pairing, there exist in these materials lower frequency 
phonons with frequency ujq which may bring us to the 
regime of intermediate ujq/A. It will be interesting to 
search for the isotope effect on the penetration depth in 
these materials (provided they are in the clean limit). 
In addition, the new class of Fe-based superconductors 
have a rather large ratio of A/lue and should exhibit the 
isotope effect on the penetration depth according to our 
theory. 



G p (cj n + uj m )G p (uj. 
1 



47TI 



C 



dz tanh — G p {u n - iz)G p (—iz), (Al) 



where the contour Co surrounds poles of tanh ^ that lie 
on the imaginary axis (see FIG. [6]). Having in mind sub- 
sequent integration over £ p , it is convenient to represent 
the electron Green's function (fTTj) as^: 



G p (iuj n ) 



2e(iuj n ) 



Q(iuj n ) + T 3 e(iuj n ) + t 1 A(iu) n ) 



e(iui n ) - £ p 



+ 



uj(iu! n ) — T 3 e(iu> n ) + t A(iu n ) 



e(iui n ) + £ P 



where e(lu) is defined as 



A 2 (w). 



(A2) 



(A3) 



Originally G p (iuj n ) in Eq. (|A2[) is defined only at the dis- 
crete set of Matsubara frequencies on the imaginary axis. 
Analytic properties of G p (z) in the complex plane are re- 
lated to the properties of function e(z). We suppose that 
on the imaginary axis e(z) is a well defined function. Let 
us denote a point where e(z) has an essential singularity 



8 



as oja- Wc argue that ua is real. Indeed, on the real axis 
the self-energy correction has a non-zero imaginary part 
only for frequencies above the threshold, u> > A + lue- 
Therefore, wa is real and u\ = A + 0(A), where 0(A) 
- corrections of order of A. Drawing branch cuts from 
points ±wa to infinity (they are shown in red zigzags in 
FIG. [5]) wc make e(z) a well defined function in the whole 
complex plane. 

Retarded and advanced Green's functions G R ' A (z) arc 
given by the value of G p (iui n ) at the upper (lower) side of 
the branch cut. They are obtained from Eq. (|A2[) by re- 
placing functions w(z), A(z), and e(z) with their analytic 



continuation to the upper (lower) side of cuts. These are 
denoted as Cj R ' A (z), A R ' A (z), and e R ' A (z) correspond- 
ingly, where first two functions were calculated in the 
Section Hi A[ and last function is defined as 

e a (z) = sign z ^J(Cj a (z)) 2 - (A«(z)) 2 , ( A4) 

where a = R, A. Now we deform the contour from Co to 
C\ (see FIG.O and analytically continue to the real ex- 
ternal frequency u> n — > —iuj. As a result, the Matsubara 
sum (|A1[) is written as 



1 poo 

^G P K + W ,„)G P K) = - J dzt a nh—[(G R (z)-G A (z))(G R (z + u;)+G A (z-oj))+ 

(G*(-z) G A (-z))(G R (-z + w) + G A (~z w ))}. (A5) 
The trace of the product of two Green 's functions G a (z 1 )G {z 2 ) with a,(3 = R,A can be written as: 

g^(z u z 2 ) + l g^(z u z 2 ) + l 



g°*P(z 1 ,z 2 )-l , g^(z 1 ,z 2 )-l 



where g a ° are the structure factors introduced ini* 



9 H {z llZ2 ) 



o) a (zi)a;^(z2) + A°(zi)A^(z 2 ) 
e a ( Zl )eP(z 2 ) 



(A7) 



While integrating over £ p , we close the integration con- 
tour in the upper half plane and use representation (|A6[) 
along with the property lme R (z) > 0, lme A (z) < 0. 
For example, integration of the product of two retarded 
Green's functions from Eq. (|A5[) yields 



d£ptr [G R (z)G R (z + u)] =-2m 



. g [z, z + u)-l 



e R (z) +e R (z + uj)' 
(A8) 

Gathering contributions from all terms and using the 
symmetry e R (— z) = —e (z), we have: 



Q(w) 



dz 



z + ui 
tanh B 



2T 



tanh — A 
2T 



dz tanh ■ 



2T 



-B, (A9) 



where 



A 



B = — 



] RR (z + uj,z) 



g RA (z + L0,z) + l 
e R {z + lu)+ e R (z) + s R (z + u) - s A (z) ' 
g AA (z + uj, z) - 1 g RA (z + uj,z) + I 
e A {z + uj) + e A {z) + e R (z + uj)-e A (z)' 



(A10) 
(All) 



We restrict ourselves to zero temperature. B — A is real, 
thus it does not contribute to the imaginary part of the 
response kernel, and consequently to cr'(oj), 



ne 2 Q"(u) 



(A12) 



The region of integration from -wa to wa in the second 
integral in Eq. (|A9|) also gives a real contribution and can 
be omitted. Finally, neglecting the difference between u>a 
and A that gives higher order corrections in A, finally we 
have for Q"{oj): 



1 



Q"(ui) = - / dzlmB 



(A13) 



which coincide with Eq. (pRJ|) used in Section Hi Bl 

In order to calculate the conductivity we use the imag- 
inary part of Cu R ' A (z) and A R ' A (z), Eqs. dHJ)- (HTJ) . Wc 
work to leading order in the coupling constant, therefore 
we expand B in Eq. (|A13[) and extract its imaginary part. 

An imaginary contribution to B comes only from the 
imaginary corrections to the self-energy, ui "(z) and 
A R "(z), Eqs. ([2l ]) - ([2"2]> . They arc proportional to the 
coupling constant, thus we expand in them. First, we 
expand e R - A (z): 



z) = e(z) sign z ± iT R (z), 



(A14) 
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with short-hand notation Using the expansion of e R ' A (z) we expand the structure 

factors and numerators in (|A11|) . After some calculation 

e(z) = ^z 2 - A 2 , (A15) we get 

rR = _ x g M >A+^ M - M -A ! ) 

2Vz 2 - A 2 v / (\z\-lo d ) 2 - A 2 

(A16) 



1 T r, A 

— Im is = T^n 7 

2 uje(z)e(z + uj) 



ACj r "{z) + z~A R "{z) AQ R "(z + lu) + {z + u)A R "(z + u) 



e(z) 



e(z + w) 

w 2 e(z)e(z + w) L ' J v ' 



To simplify this further we use the explicit form of Q R (z), A R (z), and T R (z). Inserting ImB into Eqs. (|A12[) - ()A13[) 
and changing integration variables yields 



ne 2 ttloeX H~ a {-t + 7 )(t - 7 ) + A 2 

m ^ 7_ 7+A y/{t + 7) 2 - A 2 ^(t- 7) 2 - A 2 



(A18) 



where 7 = (u> — ljd)/2. This integral can be expressed through elliptic functions, resulting in Eq. ([33 
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